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ABSTRACT 

(N ; 

t^v We study the distribution functions of mass and circular velocity for dark 

matter halos in N-body simulations of the ACDM cosmology, addressing red- 
shift and environmental dependence. The dynamical range enables us to resolve 
' subhalos and distinguish them from "distinct" halos. The mass function is com- 

pared to analytic models, and is used to derive the more observationally relevant 
y— i ■ circular velocity function. The distribution functions in the velocity range 100- 

500 km s™ 1 are well fit by a power-law with two parameters, slope and amplitude. 
\ We present the parameter dependence on redshift and provide useful fitting for- 

mulae. The amplitudes of the mass functions decrease with z, but, contrary to 
naive expectation, the comoving density of halos of a fixed velocity ~ 200 km s -1 
actually increases out to z ~ 5. This is because high-2; halos are denser, so a 
fixed velocity corresponds to a smaller mass. The slope of the velocity function 
at z = is as steep as ~ —4, and the mass and velocity functions of distinct halos 
steepen with increasing z, while the functions of subhalos do not steepen with 
z, and become even flatter at z > 2. A simple observable prediction is that the 
slope of the velocity function of isolated galaxies is steeper than that of galaxies 
in groups by as much as unity, reflecting the density biasing of high-velocity ha- 
. los. We confirm that the Press-Schechter approximation typically overestimates 

the halo mass function by a factor of ~ 2, while modified approximations provide 
?-j ' improved predictions. 
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1. Introduction 

Understanding the evolution of structure in the universe and its relation to initial con- 
ditions is a fundamental issue in cosmology. Galaxy formation involves physical processes 
such as star formation, supernovae feedback, and dust extinction which we still are far from 
understanding in full. On the other hand, the clustering of the collisionless dark-matter 
(DM) component is well understood. DM particles cluster into halos, which subsequently 
accrete more mass either by gradual infall or by merging with other halos. An accurate de- 
scription of the distribution of DM halo mass (the mass function) and its temporal evolution 
are crucial for our understanding of the evolution of the luminous galaxies that reside in 
these halos. 

A popular approximation to the mass function is provided by the Press-Schechter for- 
malism (Press & Schechter 1974, hereafter PS). The initial density fluctuation field, smoothed 
on some comoving scale R, is assumed to evolve according to linear theory until it reaches 
a critical value, 5 C = 1.69, at which time it is assumed to have collapsed and virialized. 
The critical value is determined by a linear extrapolation to the collapse time as set by the 
spherical collapse model. Many of the simplifying assumptions behind the PS formalism do 
not seem to be physically motivated, yet it turns out to predict with reasonable accuracy the 
mass function as seen in simulations; the neglected effects tend to cancel each other (Monaco 

1998) . 

While being very useful for qualitative and semi- quantitative analyses, the accuracy of 
the PS approximation may not be enough when quantitative details are concerned. For exam- 
ple, it has been found based on simulations that for halo masses below M = 10 13 — 1O 14 /i _1 M 
(depending on the specific cosmological model) the PS approximation overestimates the mass 
function by as much as a factor of ~ 2 (Gross et al. 1998; Lee & Shandarin 1999; Somerville 
et al. 2000). At the high-mass end, M £ lO 15 /i _1 M , the PS formalism underestimates the 
abundance of DM halos, especially at high redshift (e.g., Governato et al. 1999, Somerville 
et al. 2000). Improved approximations have been proposed in order to correct these inac- 
curacies, e.g., a modification of PS by Sheth and Tormen (1999, hereafter ST), practically 
replacing the spherical collapse model with an ellipsoidal collapse model (Sheth, Mo, & Tor- 
men 2000). Lee & Shandarin (1998, hereafter LS) also use a non-spherical approach, based 
on the Zeldovich approximation (Zeldovich 1970), to improve the PS approximation. 

The need for a more accurate description of the mass function is further highlighted by 
the development of a useful approach for studying galaxy formation and evolution based on 
semi-analytical models (Kauffmann et al. 1999; Somerville & Primack 1999; Baugh et al. 

1999) . Semi-analytical models try to circumvent the complications associated with baryonic 
processes by simplified prescriptions and artificially embedding galaxies within DM halos, 
which allows them to subsequently predict observable statistical quantities such as the galaxy 
luminosity function or the Tully-Fisher relation. Some of the important processes governing 
the evolution of galaxies in semi-analytical models are directly related to the DM clustering 
properties; for example, the mass function and merger histories of the halo populations are 
key ingredients in semi-analytical models (Kauffmann et al. 1999; Somerville & Kolatt 1999). 
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These need to be followed with a better accuracy than provided by the PS approximation. 
Furthermore, the PS formalism does not address substructure within halos. When two halos 
merge, the PS approximation immediately labels them as one virialized halo. However, high- 
resolution simulations (e.g., Klypin et al. 1999a, Ghigna et al. 1998) show that substructure 
does maintain its identity after merging into larger halos. This means that galactic halos 
are expected to survive tidal stripping for some time, with possible important consequences 
for collisions and starbursts (Kolatt et al. 1999, 2000). 

An alternative to modeling the complicated physics needed to predict galaxy luminosi- 
ties in individual halos is to calculate statistical properties such as the distribution function of 
the circular velocity of halos — the velocity function. By using observed luminosity-velocity 
relations, one can then relate the predicted velocity function to observational luminosity 
functions (Gonzales et al. 2000; Bullock et al. 2000b) and address halo number counts di- 
rectly. In addition, because modeling luminosities of the high-redshift galaxies is even more 
uncertain than modeling the local population, the redshift evolution of the halo velocity 
function provides a much-needed, and more direct handle on the galaxy population than the 
corresponding mass function. 

In order to obtain the velocity function we use a high-resolution simulation which has the 
force resolution and mass resolution necessary to determine the maximum circular velocities 
of halos, and is capable of resolving substructure within halos (Kravtsov, Klypin & Khokhlov 
1997). The simulation is complemented by a halo finding algorithm which also classifies halos 
into nesting levels of halos within halos (Bullock et al. 2000a; Bullock 1999; see §0). Our halo 
finder /classifier (HFC) fits each halo by an NFW density profile (Navarro, Frenk & White 
1996) and automatically assigns a circular velocity profile to each halo. For the purpose of 
comparing to observations, some of the halos are assumed to be galactic halos based on a 
simple prescription following basic observational constraints. 

In §0we present our method for extracting the mass (§ |2.1p and velocity (§ j2.2| ) functions 
for the different populations of halos (e.g., subhalos that reside in massive hosts, halos of 
galaxies in groups, etc.), while correcting for incompleteness in the halo-finding algorithm by 



an iterative procedure. In §[0] we present a simple scheme for identifying "galactic" halos. 
In §0 we compare results for distinct halos to the PS, ST, and LS predictions, and extend 
our investigation to the mass and velocity functions of substructure as well. The results for 



halos and subhalos at z = are discussed in §3.1, and their redshift evolution is presented 



in §3.2. Results for galactic halos in different environments are outlined in §3.3. We discuss 



our results and compare them to other studies in 



2. Method 

Advances in numerical simulations of the collisionless DM component allow a dynamical 
range large enough for studying substructure within halos in a volume that approaches a 
fair cosmological sample of the galaxy distribution (Klypin et al. 1999a). Using an adaptive 
refinement tree code (Kravtsov, Klypin & Khokhlov 1997), we have simulated the currently 
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popular ACDM model within a comoving periodic box of 60 h Mpc, with a force resolution 
of /res ~ 2h~ 1 kpc in the dense regions. The cosmological model is spatially flat, with 
matter and cosmological constant contributions of Q m = 0.3 and Q\ — 0.7 at z = 0. The 
Hubble constant is h — 0.7 (Hq = 100/ikms" 1 Mpc -1 ), and the fluctuation amplitude 
today is normalized by og = 1.0. The 256 3 simulated particles imply a mass resolution of 
m p = 1.1 x IO^Mq. 

Within this simulation it is possible to resolve halos inside halos to four levels down 
the hierarchy. This is done via a halo finder/classifier, based on the bound density maxima 
method (Klypin & Holtzman 1997), which has been specifically developed for the purpose of 
analyzing substructure. A detailed description of the HFC can be found in Bullock (1999) 
and in the appendix of Bullock et al. (2000a). Below are some of its basic ingredients that 
are relevant to the mass and velocity distribution functions. 

A key feature of the HFC is that it models the radial density profile of each halo with 
a universal functional form. We use here the NFW profile: 

P(r)= (r/R s )(l + r/R s y ' (1) 

where the two free parameters are an inner density p s and a characteristic scale length R s . 
It is also convenient to define the virial radius of a halo by 

MSP 



vir 



where M V - 1T is the virial mass, and virialization is defined by the mean density inside R 
being A v ; r times the mean density p. The value of the virial overdensity A V i r comes from the 
top-hat collapse model; it is about 200 for an Einstein-deSitter cosmology, and A vir ~ 340 
at z=0 for the ACDM model simulated here. The corresponding virial velocity is defined 
by Vtfj = GM vir / R vir . The profile fits were performed in the range r > 0.02i? vir . The 
fitting procedure automatically provides errors for the parameters p s and R s , which we then 
translate to errors in the other useful halo parameters such as M vir and V v i r . 

The modeling of each halo by a smooth spherically-symmetric profile helps the iden- 
tification of halos within halos. Halos at the top of the hierarchy (not within other halos) 
are termed "distinct" , while all other halos are termed "subhalos" . We only try to fit halos 
with modeled mass equivalent to at least 50 particles. As will be seen below (§ |2.1|) , the HFC 
finds halos with M ^ 2 x 10 /t M at almost 100% efficiency, while at lower masses the 
constructed halo catalogs gradually become incomplete. At z—0, a total of ~8000 DM halos 
are identified, 90% of which are distinct. Their numbers drop by approximately an order of 
magnitude for each level down the hierarchy. 

In the following two subsections we present our method for reconstructing the mass and 
velocity functions of any subset of the halos identified by the HFC. An iterative procedure 



corrects for the incompleteness at small masses ( §2.1| ). The velocity function is recovered in 



272} based on the derived mass function and the NFW profile fit. 



- 5 - 



2.1. 



Derivation of the Mass Function 



A simple count in bins of the halo population is prone to two main sources of systematic 
error. One is due to the incomplete efficiency of the HFC at masses below ~ 2 x 10 11 h^M®, 
which causes an undercount in the low-mass bins. The other is Malmquist bias due the 
low-mass cutoff imposed at 50 particles and the gradient in the mass function. In order 
to obtain an unbiased mass function, we pursue an iterative procedure as described in the 
following eight steps: 

[1] The mass M is the modeled virial halo mass M vir . We first impose a mass cutoff 
corresponding to 50 particles, i.e., M min ~ 5.5 x I0 10 h~ 1 M Q , and then count the halos in bins 
of constant logM width. The logarithmic bin size, which is about 0.4 but slightly varying 
from case to case, is determined such that we have about 15 bins in the available mass range. 
We denote these raw counts by Nj. They are shown in Figure [I] for the case of distinct halos 



[2] We evaluate the errors in each bin, ANj, as the sum in quadrature of the Poisson 
error due to the finite number of halos in the bin and the error in the assignment of halos 
to bins because of the uncertainty in the halo mass as determined by the HFC. In order to 
evaluate the latter, we produce 20 synthetic halo catalogs in which we perturb each mass 
of the simulated sample by an amount drawn at random from a Gaussian distribution of 
width equal to the corresponding HFC error. We count the halos in bins for each perturbed 
catalog, and take the error corresponding to each bin to be the standard deviation over the 
20 perturbed catalogs. The total errors, ANj, are shown as error bars in Figure [1]. 

[3] We compute the completeness function of the HFC. It is assumed to depend only 
on the number of particles in the halo, namely M, and not on redshift, or on whether 
the halo is a distinct halo or a subhalo. This assumption (which is validated a posteriori) 
allows us to evaluate the completeness only once, e.g., for distinct halos at z=0, and use 
it as is for all others cases. We find that in an extended range below 10 14 /i _1 M o , the 
mass function resembles a power law, while it drops sharply only near the low-mass end 
where incompleteness effects are manifested. We therefore assume that the underlying mass 
function is a power law all the way down to M min and determine this power law by a fit that 
uses the bins and errors determined above but excludes the first two bins near the minimum 
mass cutoff. The completeness in the j-th bin, Cj, is defined by the ratio of the raw count 
Nj in that bin and the count corresponding to the power- law fit. The completeness function 
estimated this way is shown in the inset of Figure |l|. 

[4] The binned data corrected for incompleteness, Nj/Cj, are fitted by a (temporary) 
Schechter function, according to which the number density of halos in the mass range 
(M, M + dM) is given by 



at z — 0. 



dn = $(M) dM = $*M Q e 



— M 



dM , M = M/M< 



(3) 



or 



dn = 0(M) d(logM) = 0*MV A/ d(logM) , 



(4) 
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where a = a + 1 and 0* = In 10. The fit parameters are a, 0*, and M*, and the minimum- 
X 2 fit procedure also provides the error in each one. This completes the preliminary stage of 
the analysis. Next, we correct the Schechter function for biases via an iterative procedure. 

[5] Given an assumed "true" mass function in the n-th iteration n (derived as above 
for the first iteration), we produce 20 mock catalogs of halo masses drawn at random from 
this Schechter distribution. We perturb the masses by a random deviant that is drawn from 
a Gaussian distribution of errors appropriate for that mass range, which we have derived 
before the iterations from the actual error distribution within the corresponding mass bin. 
We then bin each perturbed mock dataset and fit it by a Schechter function. A "biased" 
Schechter function 0b is obtained by averaging the Schechter parameters over the 20 noisy 
mock catalogs (and the errors in the Schechter parameters are determined by the correspond- 
ing standard deviations over the mock catalogs). Thus, the Malmquist biases lead from a 
"true" 0„ to an "observed" 0b- 

[6] To check how close we are to the desired solution, we pursue a comparison of the 
current model and data in the "observational" plane. For this, we multiply the biased 
Schechter counts by the completeness factors Cj, and obtain biased and incomplete model 
counts in bins, Nj. The measure used to evaluate convergence to the raw data is the weighted 
sum of residuals: S = — Nj) 2 / (ANj) 2 . Once we succeed in bringing S to about unity, 

4> n can be considered to be a good estimate of the unbiased mass function. 

[7] If S is smaller than in the previous iteration, then we correct our previous guess for 
the bias by: 

n+ l = (0„/0b)0n , (5) 

and go back to step [5] for the next iteration. 

[8] When S stops decreasing, we stop the iterations and adopt the current n as our 
unbiased mass function. If S is of order unity, this concludes our procedure. 

We comment that the completeness function as estimated in step [3] is in very good 
agreement with the completeness function as estimated independently by Bullock (1999), 
based on a comparison of our halo population with the halos found by a straightforward 
halo-finding technique that does not model the halo profiles. This agreement provides strong 
evidence that the estimate of the completeness function is correct. 

The use of the Schechter function for the mass function is motivated by the following two 
arguments. First, luminosity functions are known to be well-fit by the Schechter function, 
so if the mass-luminosity relation is roughly a power law, then the mass function is expected 
to be of a similar general form, with a power-law regime and a sharp drop at the large-mass 
end. Second, a similar function is predicted by the PS approximation when the standard 
deviation of the density field a(M) is assumed to be roughly a power law. We will see that 
our data is capable of determining only the two parameters of the power-law part, while M* 
is only weakly constrained, and we don't even try to fine-tune the fit with any additional 
parameters. 

Figure [I] shows the final, unbiased Schechter function (solid curve). Also shown are the 
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Fig. 1. — The mass function for distinct halos at z=0. The raw counts Nj and their errors are 
shown by open circles with error bars. The unbiased Schechter function is marked by the solid curve, 
and the associated biased and incomplete counts Nj are shown as open pentagons; their proximity 
to the raw counts is a measure of our success in recovering an unbiased mass function. Also shown 
are the predictions of the Press-Schechter (dotted) and Sheth- Tormen (dashed) approximations. 
Inset: The completeness function at the low-mass end. 



corresponding biased and incomplete counts, Nj (open pentagons). The good agreement 
between these counts and the raw counts Nj is an indication for the success of our procedure 
in obtaining an unbiased result. The final value of S in this case is 1.25, which is indeed 
on the order of unity. The PS approximation overestimates the simulated mass function by 
a factor of ~2 for M & 5 x 10 13 /i -1 M Q , then seems to cross over and underestimates the 
mass function for M > 1O 14 /i -1 M . The ST approximation agrees with the simulated mass 
function to within 10% for M^5 x 10 13 /i _1 Mq, and seems to underestimate the simulated 
mass function for M > 10 lA h~ 1 M Q . The simulated mass function at M > 1O 14 /i _1 M carries 
a large error due to the small number of halos in this mass range, and therefore the apparent 
discrepancy of the approximations in this range is of low statistical significance. 

Figure |2] compares the mass function for distinct halos and subhalos at z = 0.03 (which 
is the same as z = for all practical purposes). The raw HFC data are also shown, with the 
apparent completeness turn-off at the two lowest mass bins. The PS and ST approximations 
should be compared to the distinct halos' mass function. The success of the approximations 
is very similar to the z=0 case. 

The exponential turnoff set by M* is not constrained very well. In fact, the error in M* 
is typically larger than M* itself. This is clearly demonstrated in the subhalo mass function 
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Fig. 2. — Mass function for distinct halos (top) and for subhalos (bottom) at z = 0.03. Raw counts 
are marked by symbols with error bars. The curves are the unbiased Schechter-function fits. The 
Press-Schechter (dotted) and Sheth-Tormen (dashed) predictions for distinct halos are shown. 

of Figure 0, where the last bin actually shows an apparent increase rather than a drop. The 
biased model counts are not shown in this figure; their consistency with the raw counts is 
similar to what we had in Figure [|, with a similar value of S. 



2.2. Derivation of the Velocity Function 

The HFC also provides for each halo an NFW maximum circular velocity, V m3iX , and one 
could straightforwardly compute a raw velocity function. However, it would be impossible 
to correct this velocity function for biases following a similar procedure to the one applied 
to the mass function, because the incompleteness of the HFC depends on mass, and there 
is no one to one correspondence between velocity and mass. We therefore adopt a different 
approach, where the biases are corrected at the mass-function level, and the velocity function 
is derived from the corrected mass function. 

The conditional distribution of V max given M is not approximated by a simple func- 
tion and it strongly depends of redshift, so we choose to go from mass to velocity via the 
NFW "concentration" parameter c = R vw /R s . We draw halo masses at random from the 
corrected Schechter mass distribution, and assign to each of them a value of c drawn from 
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the conditional lognormal distribution 

~-(logc- (logc|M)) s 



P(c\M) = J- exp 



2A 2 



(6) 



where A is the conditional standard deviation of log c at a given M. The values of A and 
(logc|M) are adopted from Bullock et al. (2000a); they are different at different redshifts 
as well as for distinct halos and subhalos. The lognormal distribution assumed in @ will 
be justified a posteriori. Given M and c, the NFW profile is uniquely determined, and in 
particular 

V max = 0.465 Kir(M)[c- 1 ln(l + c)-(l + c)- 1 ] _1/2 . (7) 
We then count the halos in bins of constant log Vn iax width, and fit a power law. 

One main uncertainty in comparing the raw halo velocity function to that of galaxies is 
the response of the halos to the dissipative contraction of gas inside them. The typical effect 
is of an increase in halo maximum circular velocity, both due to the direct gravitational force 
exerted by the disk and the contraction of the halo in response to the additional inwards 
pull by the formed disk. We estimate the change in the relation between M and V m&K using 
the fitting formula of Mo, Mao, & White (1997, hereafter MMW), under the assumptions 
that (|B|) is valid, all halos have a spin parameter A = 0.05, the specific angular momentum 
of the disk is equal to that of the dark matter, and the halo mass fraction which ends up in 
the disk is ma = 0.04. The counts of the corrected V max values are also well fit by a power 
law. 

Figure [3] shows the velocity functions for distinct halos and for subhalos at z=0. Shown 
are the incomplete, raw counts as well as the corrected counts as derived from the unbiased 
Schechter mass function. The associated power-law fits are shown both before and after the 
correction for baryonic infall. Note how well this procedure corrects for the incompleteness 
at the low V ma , x end. This result is still based on the same completeness function derived in 
the previous section as a function of mass, and it thus confirms its validity as well as the 
validity of @. 

The velocity functions are clearly well fit by power laws across the whole available 
velocity range. We note, however, that a more realistic, semi-analytic treatment of disk 
masses may modify the power-law shape at the low and high velocity ends (Gonzalez et 
al. 2000). 

Also shown in Figure |3] is a simplistic, PS-based theoretical prediction for the velocity 
function, not corrected for baryonic infall. It refers to distinct halos, assuming that the 
mass function is given by the PS prediction, and that the mass and velocity are related via 
the same concentration parameter for all halos, determined by the average (logc|M). This 
approximation fails to reproduce the velocity function of the halos in our simulation by a 
factor of ~ 2 in amplitude and by a wrong slope. 
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Fig. 3. — Velocity function for distinct halos (circles, thick lines) and subhalos (triangles, thin lines) 
at z=0. Open symbols refer to the raw, incomplete counts, while filled symbols with error bars 
refer to the corrected counts derived from the mass function. The solid lines are the power-law fits 
to the corrected counts, while the dot-dashed lines are the fits to the counts corrected for baryonic 
infall. The dotted line is a prediction based on the Press-Schechter formalism, not corrected for 
baryonic infall. 

2.3. Identification of "Galactic" Halos 

In order to directly compare our simulation results with data based on observed galaxy 
rotation curves, we should identify a fraction of our HFC halos as "galactic" halos. We adopt 
the simplified assumption that every halo with mass M < lO 13 /i _1 M hosts a galaxy. The 
mass cutoff is placed to exclude groups and clusters from the galaxy count. Figure § shows 
the two-point correlation function of these "galactic" halos; it is indeed in good agreement 
with the correlation function measured for galaxies in the APM survey (Baugh & Efstathiou 
1993) for y max > 150 km s _1 . 

Further division into "galaxies" in groups and clusters versus isolated, field "galaxies" 
may yield interesting theoretical predictions that can be directly confronted with obser- 
vations. We identify isolated galaxies with the distinct "galactic" halos that do not host 
subhalos, and grouped galaxies with the "galactic" halos that do contain subhalos or that 
are subhalos themselves. Observationally, the grouped galaxies are 40 to 60% of the total 
galaxy population (e.g., Ramella, Pisani & Geller 1997; Zabludoff & Mulchaey 1998). If we 
define a halo as a subhalo only if its center lies within the i? v i r of a larger halo (as done 
originally by Kolatt et al. 2000), then the above classification scheme yields only ~ 15% in 
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Fig. 4. — Two-point correlation functions for simulated "galactic" halos with 
(symbols and dashed line) and for galaxies in the APM survey (solid line). 
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grouped "galaxies" at z = 0. When we relax the proximity requirement of subhalos and 
hosts to 3-R v i r , we obtain in the simulations a more reasonable fraction of ~ 40% in grouped 
"galaxies" . Indeed, the observer classification tends to include loose groups whose DM halos 
may not overlap. 

For the two populations of "galaxies" we assign velocities according to the mass - 
concentration relation of distinct halos. This may introduce a small error in the case of 
grouped "galaxies", which include some subhalos. 



3. Results 
3.1. Distinct Halos and Subhalos at z=0 

Figure ||] shows the difference between the mass functions of distinct halos and subhalos 
at the present epoch. | Apart from the difference in normalization, the two populations seem 
to follow a similar power-law distribution, with a slope of a = —0.83 ± 0.03 for the distinct 
halos and a = —0.82 ± 0.13 for the subhalos. As mentioned above, the turnoff scale M* 



6 We actually use z=0.03 rather that z=0 in order to compare neighboring redshift outputs (see §2.1) and 
to avoid redundancy with Figure [l]. 
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is not well constrained, with a large formal error that is mainly due to the large Poisson 
errors in the high-mass bins. The relative errors in the fit parameters grow only slightly with 
increasing z. For example, the relative error in the slope increases by ~20% between z=0 
and z=3. 

We fit the velocity function with a power law: 

dn = *(V^)dV^ = **v£ m dV max , (8) 

or 

dn = i>{V max ) d(log V max ) = d(log VW) , (9) 

with j3 — (3 + 1 and ip* = ^ In 10. We do not attempt to fit a Schechter function because 
there is not even a hint for a break at large velocities; the weak signature of an exponential 
break seen in the mass function is completely erased in the velocity function by the scatter 
in the concentration parameter (eq. |J). 

Figure |3| shows that the main qualitative results of the mass function carry over to the 
velocity function, namely, an order-of-magnitude difference in the normalization but very 
similar slopes for distinct halos (/3 = —2.82 ±0.03) and subhalos {(3 = —3.04 ±0.10). This is 
despite the fact that the corresponding distributions of concentration parameter [P(c|M)] are 
different. The fits are quite robust, with the relative errors in the slope and the normalization 
(defined at V max = 300kms _1 ) limited to a few percent, at all redshifts in the range 0-3. 

The effect of applying the simple MMW baryonic infall prescription with mj = 0.04 is 
a logarithmic shift of the halo population to larger velocities (as indicated by the arrows). 
The corresponding increase in normalization is a factor of ~ 2, both for distinct halos and 
subhalos. The change in slope, to /3 = — 2.79±0.04 for distinct halos and j3 = — 2.92±0.09 for 
subhalos, is statistically insignificant. This was not expected a priori because the baryonic- 
infall correction depends on halo concentration c (halos of higher c are assumed to accrete 
more baryonic mass) and it is therefore expected to vary as a function of mass and velocity. 
Furthermore, the correction is virtually the same for distinct halos and subhalos, despite the 
fact that the relation between mass and concentration is somewhat different for these two 
populations (Bullock et al. 2000a). A more sophisticated treatment of infall, which would 
take into account effects such as the lower fraction of baryons involved in forming the disk 
inside low-Vmax halos, may yield a change in shape (e.g., a slight turn-off at low velocities, 
see Gonzalez et al. 2000). 



3.2. Redshift Evolution of Distinct Halos and Subhalos 

The mass functions for distinct halos at three different redshifts in the range to 4 are 
shown in Figure |5|. Shown for comparison are the predictions of the PS and ST models. At 
z=0, we recover the already known result (e.g., Gross 1997, Somerville et al. 2000) that at 
low masses the PS approximation overestimates the mass function by a factor of ~2, while 
at high masses it is an underestimate. The ST approximation fares better in general, and 
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Fig. 5. — The fitted mass function for distinct halos in the simulations (solid) at three different 
redshifts (z=0, 2.3, 4, from top to bottom). The fits are truncated to span a mass range defined 
by bins for which Nj/ANj > 2. The approximations of PS (dotted) and ST (dashed) are shown 
for comparison. 

especially at low masses; it is accurate to < 10% up to M ~ 5 x 10 13 /i _1 M o at z — 0, and 
M ~ 2 x 10 12 h~ 1 M Q at z = 4. The crossover mass, where PS agrees with the simulated 
mass function, shifts from ~ 1O 14 /i _1 M at z=0 to ~ 10 12 h~ l M & at z=4. By comparison, 
the typical collapsing mass defined by er[.M*(,z)] = 1.69 (where a[M] is the linear rms 
density fluctuation on a scale corresponding to mass M), is equal to ~ lO^/i^M© at z=0 
and to ~ 1O 9 /i -1 M at z=4. 

The redshift evolution of the mass functions for both distinct halos and subhalos are 
displayed in Figure ||]. We show the evolution of the Schechter fit parameters a and 4>(M) 
at M — 10 12 h~ 1 M Q . We present the normalization by the latter rather than by 0* because 
the large scatter in M* induces a non-negligible scatter in 0* (M* is determined by the value 
of (f> at M*), and because at this relatively low mass we have good statistics throughout the 
studied redshift range. The values of M* and 0* are given in Table 1. 

The top panels of Figure ^ show the evolution in the slope at of the mass function of 
distinct halos and subhalos. A marked difference is seen between the two kinds of halos. 
While the distinct-halo mass function steepens with redshift from at = —0.85 ± 0.04 at z=0 
to a — —1.3 ± 0.1 at z — 5, the subhalo slope has a broad minimum of a = —0.9 ± 0.15 
at z ~ 1 and it then rises to a=— 0.5 ± 0.15 at z—5. A useful functional fit that traces this 
evolution is provided by the quadratic polynomial 

ac(z) = a + bz + cz 2 , z = log(l + z) . (10) 

these fits are shown in the figure, and the corresponding parameters a, b, and c are listed in 
each panel. 
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Fig. 6. — Redshift evolution of the mass functions of distinct halos and subhalos. The top panels 
show the evolution in the slope of the mass function (left - distinct halos, circles; right - subhalos, 
triangles), and the bottom panels show the overall normalization at M = 10 12 Ii~ 1 Mq. Functional 
fits to the redshift evolution are shown as solid lines, where a, b, and c are the fit parameters (see 
text). Also shown are the Press-Schechter (dotted), Sheth-Tormen (dashed) and Lee-Shandarin 
(long-dashed) predictions (relevant only for distinct halos). 

The evolution of the normalization of <fi at 10 12 /i _1 Mpc is shown in the bottom panels 
of Figure |6|. Despite the big difference in amplitude, the trends with redshift are similar for 
distinct halos and subhalos, both remaining roughly constant between z = and ~ 1, and 
then gradually declining, by roughly an order of magnitude at z=5. Indeed, the functional 
fits to the two kinds of halos show a difference of ~ 1 in the offset a and only 15% in the 
curvature c. The error in the normalization at 10 12 h~ 1 Mpc is much smaller than the error 
in 0*. The error is the sum in quadrature of the random error associated with the count in 
the bin containing M = 10 12 /i _1 Mpc and the systematic error estimated by the change in 
the value of the fitted function across the bin. 

For distinct halos we show in comparison the evolution of the PS model predictions, 
fitted at each redshift by a Schechter mass function in the range 10 11 h~ 1 M & < M < 
1O 15 ^ _1 M . The PS approximation overestimates the normalization at 1O 12 /i _1 M for z < 2 
by a factor of ~ 1.7, while at z=5 it is an underestimate by a factor of ~ 2. The slope a 
is well approximated up to z ~ 2.5, but is underestimated at higher z by as much as unity. 
The ST and LS approximations fare better than PS. They are quite similar, predicting ac- 
curately the slope up to z ~ 3. At 3 < z < 5 both approximations predict a slope steeper by 
~20% than what we find in the simulations. The ST and LS models also improve on PS by 
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Fig. 7. — The conditional mass function: conditional probability density P(M su b| Mh os t) a t subhalo 
mass of M su b = 10 12 /i _1 Mq. The host masses are divided into three mass bins: Mh os t < 1O 13 /i _1 M 
(small triangles) 1O 13 /i _1 M < M host < 1O 14 /i _1 M (medium triangles), and M host > 1O 14 /i _1 M 
(big triangles). The solid lines are linear fits to the data. 



predicting the normalization (at 1O 12 /i _1 M ) more accurately, to within 25% in the range 
< z < 3. The discrepancies in the predicted slope (at high z) and normalization (at low 
z) are small relative to the PS model, but are, nevertheless, statistically significant. 

We are also able to investigate the properties of a conditional mass function,[] denned as 
the probability density P(M su b|Mh os t) of finding a subhalo of mass M su b inside a host halo 
of mass Mhost- The conditional mass function was derived along the lines of the prescription 
described in § |2.1| . The conditional mass functions were fit to a Schechter function at each z 
output, and the fit parameters are given in Table 2. We have divided the host masses into 
three mass ranges bordered by 10 13 and 1O 14 /i _1 M , and compared the conditional mass 
functions in these three ranges. The slopes of these functions are found to fluctuate in the 
range —0.9 < a < —0.6, with no obvious correlation with the host mass. There is also no 
evidence for redshift evolution of the slopes (Table 2, second column), but note the large 
Poisson errors. 

Figure [7| shows the conditional probability density P(M sub |M host ) at M sub = 1O 12 /i _1 M , 
for the three ranges of Mh ost , as a function of redshift. The normalization here is by the 



7 Not to be confused with the temporal conditional mass function, as in the extended PS formalism, which 
is the probability of a progenitor of mass Mi at z\ given a halo of mass M2 at Z2 (e.g., Lacey & Cole 1993). 
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Fig. 8. — The redshift evolution of the velocity function for distinct halos (left panels) and subhalos 
(right panels): The top two panels show the evolution in the slope of the velocity function, and the 
bottom panels show the overall normalization at V maiX = 300 km s" 1 . Functional fits to the redshift 
evolution are shown as solid lines, and a, b, and c are the fit parameters (see text). 



number of host halos (rather than by volume), so that an integration over all subhalo masses 
per given host-halo mass yields a probability of unity. The linear fits to the data exhibit very 
little evolution, in concordance with the results for the whole subhalo population. We see no 
evidence for M host dependence in the evolution of the subhalo mass function. This could be 
useful for constructing synthetic subhalo populations for semi-analytical simulations. Since 
host mass correlates strongly with the subhalo environment density (e.g., top-hat smoothed 
at 1.5/i Mpc, Bullock et al. 2000a), Figure [?| also reflects the environment dependence of 
the subhalo mass function. 

Figure |8] shows the evolution of the parameters characterizing the velocity functions for 
distinct halos and subhalos. Here, as for the mass function, the normalization is characterized 
by the value of the_velocity function at a fixed velocity, of 300 km s" 1 , rather than by ip*- 
The fit parameters j3 and of equation (PJ) are listed in Table 3. We also provide fits to the 
evolution by the same quadratic polynomial as for the mass functions, equation fllQl). The 
errors are estimated as for the mass function; they are larger for subhalos because there are 
fewer of them. 

The evolution of the slopes of the velocity functions is qualitatively similar to that of 
the mass function. This is not surprising since the M — V max relation at all redshifts is close 
to a power law, M oc V^ ax , implying f3(z) ~ s(z)a(z). For distinct halos, s = 3.4 at low z 
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and it approaches s = 3 at high z as typical concentration values fall. For subhalos, s = 3.9 
at z = and it again approaches s = 3 at high z (Bullock et al. 2000a). 

On the other hand, the normalizations of the velocity functions remain roughly constant 
with redshift, unlike the decrease shown by the mass functions between z = 1 and 5. This 
interesting behavior is a combination of two effects. First, ignoring the redshift dependence 
of A vir , the typical halo density increases with redshift roughly as (1 + z) 3 , implying that 
halos of a given velocity correspond to objects of smaller masses at high z: M(V maK = 
const) oc (1 + z)~ 3 / 2 . The increasing halo density with redshift is counteracted somewhat 
by the tendency of high- z halos to be less concentrated (Bullock et al. 2000a). The net 
result is a very weak evolution in the comoving number density of halos of a given, high 
velocity. Since the slope of the velocity function becomes steeper at high z, the fact that 
the normalization at a fixed V mSuX = 300 km s" 1 is roughly independent of redshift implies 
that low velocity (V max ^ 200 km s -1 ) halos, and the galaxies that reside within them, are 
predicted to be more abundant at high redshift. 



3.3. Isolated versus Grouped Galaxies 



In this subsection we aim at a prediction that can be directly tested observationally. As 
described in § |2.3| , we crudely assume that each halo of mass smaller than 1O 13 /i _1 M has a 
luminous galaxy in it; this includes both distinct halos and subhalos. We then divide the 
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Fig. 9. — The velocity function for isolated galaxies (squares, short-dashed fits) and grouped 
galaxies (stars, long-dashed fits) at z=0 and 1.35. also shown for reference are the velocity functions 
of distinct halos and subhalos (solid and dotted lines). 
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galaxy population into two classes, of isolated galaxies and grouped galaxies, containing at 
z = about 60% and 40% of the galaxy population respectively. 

Figure |9] compares the velocity functions of these two classes, at z = and at z — 
1.35. The two functions display significantly different slopes, with (3 = —3.04 ± 0.10 and 
—2.11 ± 0.03 for isolated and grouped galaxies at z — 0, and with (3 = —3.44 ± 0.08 and 
—2.43 ±0.05 at z = 1.35. The ratio between the two slopes at these two redshifts is virtually 
the same, ~0.7. At z=3 (not shown), the difference between the slopes of the velocity 
functions is even larger: (5 = —3.96 ± 0.05 and —1.85 ± 0.05 respectively, a ratio ^0.45. It 
is worth mentioning that the fraction of grouped halos (according to the 3i? v i r criterion of 
273) decreases from 40% at z = and 1.35 to 30% at z = 3. 



In general, the velocity functions are well fit by power laws. The exception is the velocity 
function of grouped galaxies at £=1.35, which drops more steeply at V mSuX > 300 kms -1 . This 
is an artifact of our imposed strict upper bound on galactic-halo masses at HF/i^M©, which 
translates into a smoother drop in velocity due to the scatter in the mass-velocity relation 
(§ [2.2|) . The corresponding drop at z = is pushed to higher velocities. In the real universe 
we may expect galactic halos more massive than 1O 13 /i _1 M , and therefore the drop may be 
unphysical and should be ignored in the power-law fit. 

The predicted difference between the slopes of the velocity functions of isolated and 
grouped galactic halos reflects the "biasing" tendency for large halos to be more clustered; 
high-mass, high- velocity halos are more likely to have companions and to exist in groups than 
their low-mass, low-velocity counterparts. This relation carries over to luminous galaxies, if 
we assume that every halo that obeys our simple mass criterion hosts a galaxy. The predicted 
slope difference between groups and the field should be observable, both at low and high 
redshifts. 



4. Discussion and Conclusions 

We have studied the distributions of mass and maximum circular velocity of DM halos 
in a cosmological simulation of the ACDM model. The mass function is useful for semi- 
analytic modeling of galaxy formation, while the velocity function can be confronted with 
observations once a relation between halo velocity and disk rotation velocity is assumed. 
The high resolution allows us to address distinct halos as well as subhalos, at a range of 
redshifts, and to distinguish between field and grouped galactic halos. 

We find that the Schechter function, and in particular its power-law behavior over a wide 
range below the characteristic mass M*, is a good fit to the different halo distributions of 
mass and velocity, over the whole range studied in this simulation, and at all times monitored. 
This is once the incompleteness of the halo finder at small halos and the associated errors 
are properly accounted for. Each of these distribution functions is therefore characterized 
by a slope and an amplitude. The characteristic large scale where the function bends is 
not constrained properly because we do not sample enough large halos. At the low end, 
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we do not attempt to recover the velocity function below 100 km s -1 , and therefore cannot 
address possible deviations from a power-law at low velocities (Klypin et al. 1999b; Moore 
et al. 1999). At z = 0, for distinct halos, we find mass and velocity function slopes of 
a = —1.85 ± 0.4 and /3 = —3.82 ± 0.03 in equation (|3|) and equation (^) respectively. 

We measure the time evolution of the mass and velocity functions via the redshift 
dependence of the slope and amplitude parameters. The z dependences are described for 
convenience by simple functional fits out to z = 5 [Figures ^| and to be used in semi- 
analytic models of galaxy formation. The slope of the velocity function for distinct halos 
steepens with redshift like (3 = —3.9 — 2.0z + lAz 2 ± 0.2, where z = log(l + z). A particular 
prediction that should be tested observationally is that the amplitude of the velocity function 
at fixed Vmax = 300 km s -1 hardly varies with redshift, for all types of halos. This implies 
that the number density of halos with V max ~ 200 km s -1 and below is actually predicted to 
increase slightly with increasing redshift, as opposed to the naive expectation based on the 
way the mass function evolves. 

This prediction for the weak evolution of the velocity function can also be tested indi- 
rectly. For example, by combining this result with the observed evolution of the luminosity 
function of spiral galaxies, one can deduce the way the Tully-Fisher relation might evolve 
with redshift (cf. Gonzalez et al. 2000; Bullock et al. 2000b). To make this connection 
between the luminosity function and the velocity function using the Tully-Fisher relation, it 
will also be necessary to take into account the effects of baryonic infall. But this is impor- 
tant, since an inconsistency between this prediction and direct observations of Tully-Fisher 
evolution may question the validity of the ACDM model. 

Most earlier cosmological simulations, and analytic approximations such as Press - 
Schechter, overlook substructure within halos, which is bound to play an important role 
in galaxy formation. Our simulations and halo finder enable a study of the hierarchical halo 
population within a cosmological volume. We find that the subhalo distributions can also 
be fitted by power laws, though the errors are larger because there is an order of magnitude 
fewer subhalos than distinct halos. At redshifts of order 2 and beyond, we find that the 
subhalo population evolves differently from the distinct halo population. While the slope of 
the distinct-halo functions becomes steeper with increasing redshift, the slope of the subhalo 
becomes flatter at higher redshifts. On the other hand, the amplitudes of the subhalo and 
distinct-halo functions decrease with increasing redshift in a similar way. 

Several other simulations were used to study the evolution of substructure within indi- 
vidual clusters of galaxies (e.g., Ghigna et al. 1998, van den Bosch et al. 1999, Sensui et al. 
2000, Okamoto & Habe 1999). The results of Okamoto & Habe (1999), in particular, are 
quite similar to what we find. They see very little evolution of the subhalo mass function in 
the range < z < 2, with a power-law slope a = a + 1 = —0.6. This similarity is despite 
the fact that they concentrate on a constrained realization of a 3-cr density peak within an 
Q = 1 SCDM model and their mass range is quite different. 

For the purpose of semi-analytic modeling of galaxy formation, we also considered the 
substructure as a function of host-halo mass. We find that except for the natural "normal- 
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ization" trend (that more massive halos harbor a larger number of subhalos, see the study of 
the multiplicity function in Kolatt et al. 2000), there is only little evolution in the subhalo 
mass function for any given host-halo masses. These conditional mass functions can serve 
to improve the semi-analytical recipes for placing galaxies in simulated DM halos. 

Our derived velocity function for distinct halos agrees with the preliminary independent 
analysis of the same simulation by Gottlober et al. (1998). We improve on the preliminary 
analysis by correcting for small-mass incompleteness, placing meaningful error bars on the 
fit parameters, following the velocity function with higher time resolution and for a longer 
history (up to -2=5), and, in particular, extending our investigation to subhalos, and studying 
environmental dependence. The derived velocity function for halos identified as galactic 
halos, which is similar to that of distinct halos, is also in good agreement with Klypin et al. 
(1999b), who identified galaxies in the same simulation, as well as in a simulation of higher 
resolution, complete down to ~ 30kms _1 . They used a different halo finder and assigned 
circular velocities in a different way. They report a velocity-function slope which corresponds 
to P — —2.75 and log'?/'* 4.84, in good agreement with our values for distinct halos. This 
agreement is encouraging evidence for the robustness of our results. 

One of our most interesting predictions to be tested observationally is that galactic halos 
in groups should have a significantly flatter velocity function than more isolated galactic 
halos. The implication is that galactic halos at the high velocity end should show a stronger 
tendency to reside in denser environments. A detailed comparison to observations must 
incorporate the relation between the halo velocity and that of the luminous galaxy, which we 
crudely assumed here to be identical. White, Tully & Davis (1988) and Mo & Lahav (1993) 
already found hints for a correlation between galaxy velocity and local galaxy density, which 
seem to be qualitatively consistent with our predicted trend. However, a proper quantitative 
comparison is yet to be done, treating effects of incompleteness as a function of velocity, and 
considering in particular the slope of the velocity function as a function of the environment. 
This comparison should be done with larger, more complete datasets, both at low and high 
redshifts. 

High resolution N-body simulations accompanied by simple schemes for galaxy identifi- 
cation can thus provide a powerful tool for investigating galaxy formation, especially those 
aspects of the problem that are not yet properly addressed by semi-analytic modeling or 
by full hydro simulations. Future simulations of larger volumes, which will include more 
massive halos of cluster size, will permit better constraints on M*. Better statistics will also 
mean a more accurate evaluation of the subhalo mass and velocity functions, especially the 
conditional mass function, which could give us a better handle on the relation between host 
halos and their internal substructure. 

As a by product, our results can serve to evaluate analytic approximations of the mass 
function in the clustering process. At low redshift, we confirm earlier findings that the 
Press-Schechter approximation overestimates the true mass function by a factor of two at 
low masses, and underestimates it at high masses. This is in good agreement with Gross 
et al. (1998) and Gross (1997), and in qualitative agreement with Somerville et al. (2000) 
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and Lee & Shandarin (1999). This agreement between the simulations, using a variety 
of halo-finding algorithms, provides an additional confirmation for the robustness of our 
HFC in finding distinct halos. The latter differ somewhat in the value of the mass where 
the PS approximation coincides with the true mass function, probably due to the different 
cosmologies they investigated (Q = 1 SCDM and rCDM). We find that the approximations 
of PS as well as ST and LS overestimate the normalization of the mass function at z < 3 
and underestimate the slope at z > 3, but the discrepancies between the simulation and 
the ST and LS predictions are much smaller than the deviations of the PS formalism. All 
three approximations predict a higher rate of evolution for the mass function than seen in 
the simulations. 
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Table 1: Parameter fits for the mass functions 







distinct halos 






subhalos 




z 


a 






a 


log 4>l 


logM* 


0.00 


-0.85 ±0.04 


-5.37 ±0.41 


15.82 ± 16.24 


-0.72 ±0.10 


-5.35 ±0.41 


15.03 ± 16.72 


0.25 


-0.86 ±0.04 


-4.90 ±0.34 


15.23 ± 16.26 


-0.80 ±0.10 


-5.52 ±0.31 


14.95 ± 16.53 


0.70 


-0.93 ±0.06 


-5.24 ±0.36 


15.37 ± 15.96 


-0.89 ±0.16 


-5.51 ±0.41 


14.64 ± 15.56 


f.35 


-1.00 ±0.05 


-5.72 ±0.32 


15.53 ± 15.82 


-0.92 ±0.15 


-5.90 ±0.56 


14.80 ± 16.70 


f.99 


-1.06 ±0.05 


-5.91 ±0.26 


15.46 ± 15.85 


-0.85 ±0.11 


-5.80 ±0.32 


14.79 ± 16.78 


2.95 


-1.16 ±0.06 


-5.77 ±0.33 


14.92 ± 15.81 


-0.94 ±0.18 


-6.08 ±0.33 


14.54 ± 17.11 


4.00 


-1.32 ±0.09 


-6.26 ±0.23 


14.68 ± 14.20 


-0.65 ±0.08 


-5.81 ±0.42 


14.88 ± 16.85 


4.92 


-1.29 ±0.09 


-6.35 ±0.28 


14.61 ± 14.85 


-0.46 ±0.15 


-5.72 ±0.46 


15.05 ± 18.51 



a units of [ft-iMpc]" 3 
b units of Ii-I-Mq 
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Table 2: Conditional mass function fits 







M host < 1013 

nost 




z 


a 




logM* 


0.00 


-0.71 ±0.07 


-5.52 ±0.37 


14.92 ± 16.08 


0.08 


-0.80 ±0.09 


-5.57±0.28 


14.67 ± 16.08 


0.25 


-0.72 ±0.07 


-5.52 ±0.39 


14.88 ± 16.35 


0.41 


-1.06 ± 0.12 


-6.15 ±0.27 


14.35 ± 15.80 


0.50 


-0.95 ±0.11 


-6.07±0.35 


14.63 ± 16.07 


0.70 


-0.87 ±0.09 


-5.92 ±0.47 


14.81 ± 15.96 


1.35 


-0.87 ±0.07 


-5.92 ±0.26 


14.82 ± 15.97 


1.99 


-0.66 ±0.09 


-5.35 ±0.39 


14.86 ± 15.71 


2.95 


-0.98 ±0.16 


-6.03 ±0.51 


14.36 ± 17.39 




10 




14 


z 


a 




logM* 


0.00 


-0.61 ±0.08 


-5.45 ±0.21 


14.99 ± 16.70 


0.08 


-0.68 ±0.09 


-5.70 ±0.55 


14.92 ± 16.90 


0.25 


-0.91 ±0.16 


-5.90 ±0.52 


14.50 ± 16.58 


0.41 


-0.69 ±0.08 


-5.47 ±0.31 


15.03 ± 16.69 


0.50 


-0.80 ±0.15 


-6.03 ±0.48 


14.76 ± 16.69 


0.70 


-0.58 ±0.07 


-5.54 ±0.40 


15.42 ± 17.18 


1.35 


-0.42 ±0.10 


-6.00 ±0.55 


16.58 ± 17.32 






MLst > io 14 




z 


a 


log 4>l 


log Ml 


0.00 


-0.65 ±0.11 


-5.49 ±0.55 


14.87 ± 17.00 


0.08 


-0.72 ±0.09 


-5.80 ±0.53 


14.96 ± 16.57 


0.25 


-0.64 ±0.12 


-5.46 ±0.63 


15.03 ± 16.91 


0.41 


-0.71 ±0.11 


-5.94 ±0.52 


14.90 ± 17.13 


0.50 


-0.59 ±0.17 


-5.57 ±0.71 


14.93 ± 18.13 


0.70 


-0.70 ±0.10 


-5.86 ±0.45 


14.84 ± 16.32 


1.35 


-0.78 ±0.10 


-5.96 ±0.34 


14.74 ± 16.38 



a units of [h-^lpc]- 3 
b units of hr x MQ 
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Table 3: Parameter fits for the velocity functions 



distinct halos subhalos 



z 


log^* 


P 


log^C 


P 


0.00 


4.93 ±0.06 


-2.82 ±0.03 


4.22 ±0.20 


-3.04 ±0.10 


0.50 


6.25 ±0.08 


-3.39 ±0.04 


6.96 ±0.39 


-4.16 ±0.20 


1.35 


5.86 ±0.09 


-3.19 ±0.04 


4.40 ± 0.38 


-3.00 ±0.18 


2.31 


6.66 ±0.11 


-3.51 ±0.05 


4.03 ±0.27 


-2.84 ±0.13 


2.95 


6.84 ±0.07 


-3.57 ±0.03 


3.02 ±0.21 


-2.66 ±0.10 


4.00 


7.54 ±0.10 


-3.96 ±0.05 


0.76 ±0.16 


-1.62 ±0.08 


4.92 


6.00 ±0.26 


-3.33 ±0.12 






a units of [/i _1 Mpc]- 


-3 







